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Abstract 

We find all scalar second order evolution equations possessing 
an sh-valued zero curvature representation that is not reducible to a 
proper subalgebra of SI2 . None of these zero-curvature representations 
admits a parameter. 

For more than twenty years, researchers are being attracted by the prob- 
lem of classification of nonlinear systems possessing a zero-curvature rep- 
resentation (ZCR). Efforts are focused on ZCR's taking values in a non- 
solvable Lie algebra g and depending on a nonremovable parameter, in 
expectation that they will be suitable for the Zakharov and Shabat [16] 
formulation of integrability (S-integrability) and hence admit soliton solu- 
tions. But even without parameter, a ZCR may be useful for construction of 
Backlund transformations, nonlocal symmetries and nonlocal conservation 
laws. The problem of finding a ZCR is then equivalent to that of computing 
finite-dimensional linear coverings in the sense of Krasil'shchik and Vinogra- 
dov [3], which are very often just realizations of the Wahlquist-Estabrook 
prolongation structures [15]. However, existing computational procedures 
are insufficient for solving general classification problems, unless in combi- 
nation with methods based on different criteria of integrability. The most 
complete lists of integrable systems obtained so far resulted from the formal 
symmetry approach [6, 7, 10]. 

In this paper we apply the method of [4, 5] and find all second order 
scalar evolution equations 

Ut ^ F{t,X,U,Ux,Uxx) (1) 

possessing an sl2-valued zero curvature representation that is irreducible in 
the sense of being not reducible to a proper subalgebra of SI2 . We arrive at 
a single previously unnoticed class of equations parametrized by a single 
function of the coordinates t, x. We also distinguish a particular subclass of 
equations that admit a single conservation law. None of the corresponding 
ZCR's admits a substantial parameter, which is in accordance with the 
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general belief that no second order scalar evolution equation is S-integrable. 
All the previously known ZCR's [1, 8] for second order scalar evolution 
equations turn out to be reducible to a solvable subalgebra and as such fall 
outside our classification. 



1 Preliminaries 

Let -E be a nonlinear partial differential equation on a number of functions in 
two independent variables t, x. Let g be a non-solvable matrix Lie algebra. 
By a Q-valued zero- curvature representation (ZCR) for E we mean a q- 
valued one-form a = Adx + B dt such that 

da = \ [a, a] (2) 

holds as a consequence of E. 

Let G be the connected and simply connected matrix Lie group as- 
sociated with 0. Then for an arbitrary G- valued function S, the form 

= dSS~^ + SaS~^ is another ZCR, which is said to be gauge equiv- 
alent to the former. Gauge equivalent ZCR's may be viewed as identical 
geometric objects (connections). A fl- valued ZCR is said to be reducible if 
it is gauge equivalent to a ZCR taking values in a proper subalgebra of fl; 
otherwise it is said to be irreducible. 

Let us proceed to a description of the general algorithm of computing 
ZCR's [4, 5] as we use it here. For simplicity we restrict ourselves to a single 
non-linear nth order evolution equation 

Ut = F{t,X,U,Ui,...,Un) (3) 

Here t, x are coordinates, u is a single field variable, and ui = u^, M2 = u^x, 
etc., represent the derivatives. Let us consider the corresponding equation 
manifold E, that is, the infinite-dimensional space R°° endowed with coor- 
dinates t, X, u and Ufe, 1 < fc. The vector fields Dx — d/dx + ui d/du+ ■ ■ ■ + 
Uk+i d/duk + ■■; Dt^ d/dt + F d/du +■■■+ D^F d/duk + ■■■ defined on 
E generate a diffiety structure in the sense of [3] and encode all essential 
geometric information about the equation. In these terms, a ZCR for cq. (3) 
is a pair of g-valued functions A,B on E satisfying eq. (2), which may be 
written as 

DtA-D^B+[A,B]=0. (4) 

Let us introduce operators Dj acting on an arbitrary g-valued function C 
on E as follows: 

D,C = D,C-[A,C], DtC = DtC-[B,C]. (5) 

Operators Dx,Dt commute whenever {A,B) is a ZCR. We also set Di = 
Dx o • • • o Dx {i times). 

By [4] for every ZCR there exists a characteristic matrix R, which 
is a fl-valued function defined on E (see also the independent work by 
Sakovich [11]). The following proposition is proved in [4, Prop. 2.7 and 
Prop. 3.9], see also [5, Prop. 2] 
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Proposition 1 (1) The characteristic matrix R for a ZCR of the evolution 
equation (3) satisfies 



s.«=i:(-s).(IC«) 



(6) 



(2) Gauge-equivalent ZCR's have conjugate characteristic matrices. 

Using Proposition 1(2), one can restrict the gauge freedom by requiring 

that the characteristic matrix R be in the Jordan normal form. To fix the 
remaining gauge freedom, due to the stabilizer S <Z G oi R, one can further 
transform the matrix A. See Section 2 for details. 

In the sequel we consider a ZCR Adx + Bdt taking values in SI2. We 
shall write the two sl2-matrices as 



Reducible ZCR's will be excluded from our classification. A subalgebra 
of SI2 to which the ZCR may be reduced is either an abelian algebra or 
the two-dimensional solvable subalgebra rcprcscntablc by lower triangular 
matrices. Obviously, a ZCR taking values in an abelian algebra is equivalent 
to one or more conservation laws. What concerns solvable algebras, the 
situation is not much different. 

Definition 2 An sl2-valued ZCR satisfying the condition 02 = 62 = is 
said to be lower triangular. 

For a lower triangular ZCR, it follows from cq. (4) that (j) = aidx + bi dt 
is a conservation law. Let h be the potential of (j); then, by the same eq. (4), 
(p' = (03 dx + 63 dt)e^^ is a conservation law nonlocal over the potential h. 
This situation will be referred to as a chain of conservation laws. Clearly, 
one can reconstruct the reducible ZCR from the corresponding chain (0, (j)'). 
In this way, reducible sl2-valued ZCR's are equivalent to certain chains of 
conservation laws. It also follows that methods to find them are to be sought 
among methods to compute nonlocal conservation laws. 

Proposition 3 Let the matrices (7) form a ZCR for for the evolution equa- 
tion (3). Suppose that 02 = 0. Then also 62 = or the ZCR is gauge 
equivalent to zero. 

Proof Let us denote by G the matrix (4) evaluated at 02 = 0. By as- 
sumption, G is zero on the equation manifold E. If 62 ^ 0, then from the 

condition = C2 = —Dxb2 + 2ai62 wc compute that ai = ^Dx^2/&2 on 
E, and then from the condition = ci = DfOi — D^bi — 0362 we compute 
that as = |-Dtx&2/&2 ~ ^-0*62 -Da; 62/^2 ~ DJ}i/b2 on E. Let us introduce 

a bmction G by the requirement 63 = —b\/b2 — Dtbi/b2 + biDib2/b\ + 



\Dub2/bl - l{Dtb2f/bl + G/&2. Then = C3 = -D^G/b2, which in the 




(7) 
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case of an evolution equation implies that G is a function of t only. Then, 
under the above substitutions for 01,03,63, the gauge matrix 



/ 6-^/2 ^ 

sends A to zero and B to 



1 
G 



which depends on t at most. The last matrix is sent to zero by gauge 
transformation with the gauge matrix composed of independent solutions 
of the equation su = Gs. 



2 The classification 

Let us consider a second order evolution equation (1) along with the SI2- 
matrices A, B satisiying eq. (4) but not reducible to a solvable subalgebra. 
We also assume that dF/duxx 7^ 0. Following [5], we consider the two cases 
distinguished by their Segre characteristics separately. 



2.1 The nilpotent case 

Under the notation (7), the Jordan form for R corresponds to ri = 0, 
?*2 = 0, r3 = 1. The normal form for A, obtained in [5], is given by the 
single requirement ai = 0. Indeed, whenever 02 7^ (otherwise the ZCR is 
either lower triangular or trivial by Proposition 3), then one can set ai = 
in a general matrix A by means of the gauge matrix 

1 

ai/a2 1 

from the stabilizer of R. 

The equation (6) then reduces to the system Tj = 0, i = 1, 2, 3, where 

Ti := 2D J a2 + £»xa2 - + 62, 

I2 ■— ^T, ^2) 

^ dF ^ dF ^ dF dF 
73 := -Dxx^ + Dx- 2- 0203 - ^ - 26i. 

dUxx OUx (JUxx OU 

Then 02 = by the second equation and 62 = by the first equation, 
whence the ZCR is lower triangular. Consequently, this case is void in our 
classification. 
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2.2 The semisimple case 



It will be convenient to change the notation for sl2-matrices to 

^ _ / ai 02 + 03 

Va2 - as -ai 

The Jordan form for R has r2 = = with r := ri arbitrary. Unlike in [5], 
we choose the normal form for A characterized by the single requirement 
as = 0. And indeed, whenever a2 + a3 ^ 0, which is irrestrictive by Proposi- 
tion 3, one can set 03 to zero by a gauge transformation from the stabilizer 
of R. The relevant gauge matrix is diagonal with the diagonal entries h and 
1/h, where h = ((02 - a3)/(a2 + 03))^/^. 

Eq. (4) and (6) then assume the form Si = = Ti, i = 1, 2, 3, with 

Si = -Dttti + D^bi + 2a2&3, 
^2 = --0/0,2 + D^b2 - 20163, 
53 = Dxbs + 20261 - 20162, 

Ti = -Dtr - —r + dJ —r) - D^J r -4- ral 



dF 

T2 = -63 + 2- 0102, 



dux 
^F_ 

dux ^ V duxx J ^ r duxx ^ du^ 



dF f OF \ D^r dF dF ^ 
T3 = -62 + 7^ — 02 - ^F)^ \ 02 - 2 02 - D,^a2. 



(8) 



If 02 — 0, then we have 62 = 63 = by the last two equations and 
the ZCR reduces to a single conservation law. Therefore, we assume that 
02 7^ in the sequel. 

Proposition 4 As solutions to eq. (8), functions r. oi, 02, 63 cannot depend 
on jet coordinates other than t, x, u, Ux, Uxx; whereas functions hi, 62 cannot 
depend on jet coordinates other than t, x, u, Ux, Uxx, Uxxx- 



Proof We may assume, without loss of generality, that the functions r, oi, 
02, 63 depend ont,x,u, . . . , Uk and the functions 61, 62 depend ont,x,u, . . . , 
Uk,Uk+i for some k > 2. We perform a downward induction, each step 
of which consists in deriving appropriate differential consequences of the 
system (8). Thus, let k > 2. Then we have 

Q_ m ^ ^ dF dr 
duk+2 duxx duk ' 
but dF/dUxx ^ 0, whence r does not depend on u^. Then, similarly, 

o db2 

duk+1 
^ dF da2 
duxx duk ' 



dS2 _ dT3 
duk+2 duk+i 

dS2 ^ dTs 
duk+2 duk+i 
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whence 62 does not depend on Uk+i and 02 does not depend on Uk- Finally, 



„ dSi dS3 dT2 , dF dai 
-2a2^ h ^ h -w— = 4a2;^ ^— , 

2a2- h h - — = 4a2- 



' duk+2 duk+i duk duk+i 
= -2a2-^ + - ^ = 2^ 

aUfc+2 dUk 

whence ai, 63 do not depend on Uk and bi does not depend on Uk+i (recall 
that a2 7^ 0). This completes the induction step. 

Under the restrictions established in Proposition 4, the determining sys- 
tem (8) becomes an overdetermined system of partial differential equations. 
As such, it can be solved routinely, but its solution is troublesome even with 
the employment of software capable of automated deriving of differential 
consequences. The reason is that the class of second order evolution equa- 
tions is invariant with respect to a large group of contact transformations 
X = x{t, X, u, Ux), u = u{t, X, u, Ux), i = i{t). Below we shall apply a series 
of suitably chosen contact transformations to achieve substantial reduction 
of the matrix A. 

Proposition 5 For every second order evolution equation (1) possessing 
an irreducible sh-valued ZCR there exists a contact transformation such 
that the transformed 02 depends on t, x, u, Ux at most. 

Proof Let functions r,ai,bi depend on t,x,u,Ux,Uxx,Uxxx as in Propo- 
sition 4. Taking successively the derivatives dSi/duxxxx, dS2/duxxxx, T2, 

dn/duxxx, aTi/aw,,,,, ds./duxxx, m/Ouxxx, d'^Si/dul,,, a^^s/a^L. 

one may check routinely that 

d'^a2 _ , da2 dai 
—^=0 and ai- 02 = 0, 

OU^x ""^xx OUxx 

are among differential consequences of the system (8). Hence, 

(a) a2 is linear in Uxx, i-e., 02 = a2i{t,x,u,Ux)uxx + a2o{t,x,u,Ux); 

(b) the ratio ai/a2 does not depend on Uxx- 

Now, if 021 = 0, then the statement is proved. Otherwise, let /i, /2 be two 
functionally independent solutions of the linear equation 

020 df df df 

^ + "^ir + ir = ^- 

021 oUx ou ox 

In particular, both fi and /2 do depend on u^- Then t = t, x = fi, u = f2 
andwj = [d f2/ dux) / {d fi/ dux) satisfy the well-known necessary conditions 



6 



of being a contact transformation: 



du 


du 


du 


dux 


dx 
dx 


du 


dx 


dx 

+ Ux^ 

du 


dux 


dx 



Under this transformation, Adx + Bdt becomes Adx + B dt with dx 
DxX dx + Dtx dt, di = dt, so that 



A = A DxX 



df2 df2 df2 

n, "T Ux „ + Uxx n 

dx du dUx 



a2 df2 
a2i dux 



A, 



where we have used eq. (9). Hence 

021 



A = 



ai/a2 1 
1 -ai/a2/' 



df2/dux 

which is independent ofuxx, hence ofuxx, by virtue of statement (b) above. 

Theorem 6 Every second order scalar evolution equation (1) possessing 
an irreducible sh-valued ZCR is transformable to an equation of the form 



dl3 



d^p 



— U Uxx + Stt-oW Ux + 4:pUx + -r-q - 4— u'' - 4 — u 



dx 



dx 



d^p , d(3 



dx"" 



dx 



dp 



dx 



(10) 



through a contact transformation. Here P is an arbitrary function of t, x 
with dp/dx ^ 0. The ZCR is then Adx + Bdt with 




B = 



f dp ^3 d^p 

-—Ux + 4 -TT-oU 

dx u dx 



4/3- 



4/3- 



2—u 
dx 



dp 
dx 



2^y 
dx 

d + ^u 
u dx^ I 



(11) 



Proof Following Proposition 5, wc; assume that the matrix A depends on 
t,x,u,Ux at most. One may chock routinely that 

5^02 . , 9^02 , 9^02 da2 _ 

, = and — — h Ux — — — = 

du^ dxdux dudux du 

are among differential consequences of the system (8). The general solution 
is a2 = dh/dx + Uxdh/du = Dxh for a suitable function h{t,x,u). If 02 
does depend on u^, then dh/du ^ 0, whence t = t, x = h, u = x is & point 
transformation. If 02 does not depend on Ux, then h does not depend on 
u, but does depend on x (otherwise 02 = 0), and t = t, x = h, u = u \s & 
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point transformation. In both cases A = AD^h = Aa2, whence 02 = 1 in 

the transformed matrix A. 

With a2 = 1, one can check routinely that dai/dux = is among 
differential consequences of system (8). If moreover dai/du = 0, then A is 
completely independent of u and its derivatives, and then so is B, whence 
the ZCR is gauge equivalent to zero. Therefore, we shall continue with 
dai/du ^ 0. Then we can apply a point transformation x = x, u = 1/ai, 
which sends ai to 1/u (this choice prevents terms quadratic in Ux from 
appearing on the right-hand side of eq. (1)). It is then a matter of routine 
to compute all possible forms of the right-hand side F of eq. (1) and also 
the corresponding matrices B. 

There seem to be no earlier appearance of the class (10) in the literature, 
let alone its 'simplest' member ut = u^Uxx + ^xux — — 4u. 

The results would be incomplete if we do not establish irreducibility 
of the ZCR (11). Since reducibility implies existence of at least one local 
conservation law, we shall start with the following result. 

Proposition 7 Within the class (10), the only equations to possess a con- 
servation law are those with 

where p, q are arbitrary functions oft such that {pq)t ^ 0. In all these cases 
the equation has a single conservation law 

pe^- + ^ / 1 (M)t(pe2- ^ ^^-2.^ 



, 1 (p«e^^ + gte-^-)(pe^-+ge-^-) l , . 
^2 (pe2--9e-2-) u ^ ' 

{pqU3p\^^ + 2pq + 3q^e-^n 
ipe'^-qe-^n' 

Proof A routine computation shows that any characteristics V' of a con- 
servation law depends on t,x,u at most and satisfies the equations 



ox ou u at ox 



The rest is easy. 

Another computation shows that for none of the equations of the 
class (12) the corresponding ZCR (11) can be reduced to the lower tri- 
angular form with multiples of (13) on the diagonal. Thus, the ZCR's (11) 

are indeed irreducible. 

Finally, a remark on equations determining pseudospherical surfaces 
(PSS equations) is due. In anticipation of finding new S-integrable nonlin- 
ear systems, a number of attempts have been made to classify equations 
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describing pseudosphcrical surfaces (PSS equations), see [14] and references 
therein. Even though being a PSS equation is equivalent to possessing an 
sl2-valued ZCR, the classification of second order scalar evolution PSS equa- 
tions as obtained by Reyes [8] (see also [9, 2]) has no intersection with ours. 
This seeming paradox is easily resolved. In fact, each of the ZCR's found 
by Reyes is reducible to the lower triangular form (the generalized Burgers 
equation) or even to a single conservation law (the other equations), which 
are disregarded in our classification. On the other side, equations (10) are 
not integrable, hence do not enter the classification of integrable equations 
by Svinolupov and Sokolov [12, 13], which was the starting point of the 
Reyes work. 
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